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OO \ ABSTRACT 

We exploit very deep mid-IR (MIR) and X-ray observations by Spitzer and Chandra 
in the GOODS North field to identify signs of hidden (cither starburst or AGN) activity 
in spheroidal galaxies between z ~ 0.3 and 1. Our reference is a complete sample of 168 
morphologically classified spheroidal (elliptical/lenticular) galaxies with zab < 22.5 
f— ^ selected from GOODS ACS imaging. Nineteen of these have 24-^tm detections in the 

GOODS catalogue, half of which have an X-ray counterpart in the 2 Ms Chandra 
catalogue (6 detected in the 2-10 keV X-rays hard band), while about 25% have 1.4 
GHz fluxes larger than 40 /zJy Traces of hidden activity in the spheroidal population 
Q ■ are also searched for in the deep X-ray images and 14 additional galaxies are detected 

|— | | in X-rays only. The nature of the observed MIR emissions is investigated by modelling 

c/3 . their Spectral Energy Distributions (SEDs) based on the available multi-wavelength 

photometry, including X-ray, UV, optical, near-IR, MIR and radio fluxes, and optical 
spectroscopy. The amount of dust derived from the IR emission observed by Spitzer 
appears in excess of that expected by mass loss from evolved stars. 

When the available independent diagnostics are compared, in general they provide 
consistent classifications about the nature of the activity in our spheroidal population. 
Given that, in principle, none of these diagnostics alone can be considered as conclu- 
sive, only trough a panchromatic comparison of them we can reach an accurate com- 
prehension of the underlying physical processes. In particular, our multi-wavelength 
analysis of the X-ray and MIR properties leads us to conclude that at least 8 of the 
19 24-/im bright sources should hide an obscured AGN, while the X-ray undetected 
sources are more likely dominated by star formation. We conclude that ~30 objects 
(~20%) of the original flux-limited sample of 168 spheroidal galaxies in the GOODS- 
North arc detected during phases of prominent activity, of both stellar and quasar 
origin. Due to the short expected lifetimes of the IR and X-ray emissions, this fraction 
might imply a significant level of activity in this class of galaxies during the relatively 
recent cosmic epochs — z ~ 0.3 to z ~ 1 - under investigation. 

Key words: galaxies: active - galaxies: elliptical and lenticular - galaxies: high- 
redshift - infrared: galaxies. 
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1 INTRODUCTION 

The origin of spheroidal (elliptical and lenticular) galaxies 
has been, and still partly remains, a controversial issue. Pub- 
lished results from high redshift galaxy surveys appear not 
unfrequently in disagreement with each other. One example 
is the apparent conflict between reports of the detection of 



massive elliptical galaxies at very high redshifts (e.g. Glaze- 
brook et al. 2004; Cimatti et al. 2004; McCarthy et al. 2004; 
Daddi et al. 2005; Saracco et al. 2005) and the indications 
for a decline in their comoving number densities at z > 1 
(e.g. Franceschini et al. 1998; Fontana et al. 2004; Faber et 
al. 2005). 

Some recent progress in this field has been achieved 
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in particular by the combined use of the IR multi- 
wavelength coverage offered by the Spitzer Space Telescope, 
the unique imaging capabilities of HST/ACS, and the enor- 
mous photon-collecting power of spectrographs on large 
ground-based telescopes. This wealth of data significantly 
refined our knowledge of the build up of galaxy populations 
and revealed evidence of "downsizing" in galaxy formation, 
i.e. star formation ending first in massive galaxies than in 
lower mass ones, as observed both at low and high redshift 
(Cowie et al. f996; Tanaka et al. 2005; van der Wei et al. 
2005; Treu et al. 2005; Juneau et al. 2005; Bundy et al. 2006; 
Borch et al. 2006; Jimenez et al. 2006). 

In particular, Treu et al. (2005) studied the evolution of 
the fundamental plane as a diagnostic for star formation and 
mass assembly history of early-type galaxies. They suggested 
that most of the stellar mass in spheroidal galaxies formed at 
z > 2, with subsequent activity continuing to lower redshifts 
(z < 1.2). The fraction of stellar mass formed at recent times 
depends strongly on galactic mass, ranging from < 1% for 
masses above 10 11 ' 5 M Q to 20%-40% below 10 11 M . 

Various approaches have been followed in order to de- 
tect recent or ongoing activity in spheroidal galaxies to 
z ~ 1. One approach was to look for close companions or 
morphological disturbances, and how they evolve as a func- 
tion of redshift (Le Fevre et al. 2000; Patton et al. 2002; 
Khochfar & Burkert 2003; Lin et al. 2004; Cassata et al. 
2005), under the assumption that the physical mechanism 
triggering the activity and responsible for the progressive 
build-up of galaxies are interactions and merging. Another 
way of attacking the problem was to consider optical colour 
information and look for blue ellipticals to find the presence 
of young stellar populations, hence recent star-formation, 
e.g. triggered by merging (Hogg et al. 2003; Bell et al. 2005; 
Cassata et al. 2006) . As an example, Treu et al. (2005) pro- 
vided independent support for recent activity in spheroidals 
up to z ~ 1 through spectroscopic ([O II] emission, H5) and 
photometric (blue cores and broadband colors) diagnostics. 

An important related question is how the observed 
galaxy build-up matches with evidences about the growth of 
nuclear super-massive black-holes (BHs) during the quasar 
and AGN phase. A close relationship of the two is suggested 
by the observed correlations of the BH mass with the host 
galaxy mass and central velocity dispersion in the local uni- 
verse (Ferrarese & Merritt 2000) and by the concomitant 
presence of star-formation and AGN activity sometimes di- 
rectly observed in forming galaxies (Alexander et al. 2005). 
Many have suggested that AGN may be the cause or at least 
related to the quenching of galaxies (see Silk & Rees 1998, 
Croton et al. 2005, Hopkins et al. 2006, among others). 

We follow in this paper an alternative line of investi- 
gation on high-z (z ~ 0.7) spheroidal galaxies by looking 
for signs of optically hidden activity in their mid-infrared 
(MIR) and hard X-ray emissions (Rigby et al. 2004, Geor- 
gantopoulos et al. 2006). Ongoing star- formation and young 
massive stars in local galaxies appear mostly embedded in 
dusty, optically thick molecular clouds. Related to this, the 
Infrared Space Observatory (ISO) has shown the presence of 
Polycyclic Aromatic Hydrocarbon molecules (PAHs) in the 
spectra of some local early-type galaxies (i.e. Madden et al. 
1999; Xilouris et al. 2004). These results have been recently 
confirmed by observations with the Spitzer Space Telescope 
(Pahre et al. 2004, Bressan et al. 2006). However, the ori- 



gin of the dust responsible for the MIR emission in early- 
type galaxies could be attributed either to merging events 
(Xilouris et al. 2004), cooling flows (Fabian et al. 1991) or 
to mass loss from late-type stars (Knapp et al. 1992). 

In addition, a significant (even major) fraction of AGN 
accretion is also expected to be optically obscured and to 
emerge in the MIR. Then an effective way of identifying 
hidden star- formation and/or AGN activity in high-redshift 
spheroidal galaxies is to exploit deep mid- and far-IR pho- 
tometry with Spitzer MIPS on sky areas where morpho- 
logical classification is possible from deep HST/ACS high- 
resolution imaging. 

Hard X-ray observations, in turn, provide an even more 
straightforward way to identify obscured AGNs, since hard 
X-rays are typically produced in copious amounts by them, 
and are able to penetrate the dust absorbing the optical, 
UV and soft X-ray light emitted by the nucleus. In addition, 
when the sensitivity of the X-ray data is high enough, they 
bear potential independent information on deeply embedded 
young stellar populations (Ranalli et al. 2003; Franceschini 
et al. 2003). Therefore, deep hard X-ray data combined with 
MIR information can provide extraordinary tools for unveil- 
ing hidden AGNs and starbursts in morphologically classi- 
fied spheroidal galaxies and to constrain the nature of these 
emissions. 

Many recent papers have studied the presence of AGN 
in high-z early-types (e.g. Georgakakis et al. 2006, Grogin 
et al. 2005, and Brand et al. 2005). Similar works at lower 
redshift are based on the Sloan Digital Sky Survey (SDSS, 
Stoughton et al. 2002) that study AGNs in high-mass early- 
type galaxies (Miller et al. 2003, Kauffmann et al. 2004). 
More recently, Capetti and Balmaverde (2005, 2006 among 
the others) explored the connection between the multiwave- 
length properties of AGNs in nearby early-type galaxies and 
the characteristics of their hosts. Also Treu et al. (2005) 
found that ~5% of their E+S0 galaxies are identified as 
AGNs, based on an X-ray luminosity above 10 42 erg/s. An 
additional ~10% of the spheroidal sample is also detected 
with X-ray luminosities between 10 40 and 10 42 erg/s, which 
are interpreted as mostly due to a low- luminosity AGN or 
to ongoing star formation. 

In this context, the GOODS-North (GOODS-N), with 
its deepest available MIR (from Spitzer) and deeper hard 
X-ray (from Chandra) data, offers a unique opportunity of 
investigating activity phenomena in spheroidal galaxies. 

We focus in this paper on the specific comparison be- 
tween the Chandra X-ray view and the Spitzer 24 am view 
of AGN and star-formation activity (Rowan- Robinson et al. 
2005, Franceschini et al. 2005), investigating in particular if 
the Chandra and Spitzer indicators agree and correlate with 
one another when they are combined. 

Based on a flux-limited optical galaxy sample (zab < 
22.5) by Bundy et al.(2005), including 168 morphologi- 
cally classified elliptical/lenticular (E/S0) galaxies in the 
GOODS-N, we have identified those with a MIR counterpart 
in the MIPS 24 ^m GOODS catalogue (S 2 4 M m > 80fi Jy). We 
have also cross-correlated the optically-selected sample with 
the 2 Msec X-ray image of Chandra reported by Alexan- 
der et al. (2003). For these galaxies we have built a multi- 
wavelength photometric catalogue including IRAC data at 
3.6, 4.5, 5.8 and 8 /jm (from the public GOODS release), 
Subaru and HST U,B,V,R,i,z and HK' data (Capak et al. 
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2004), and radio VLA 1.4GHz data (Richards et al. 2000). 
Spectroscopic redshifts are available for all the sources of 
our sample. The galaxies are found at 0.1 1zz< 1.2, with a 
median redshift < z >= 0.68. About half of the spheroidal 
galaxies with MIR emission present an X-ray counterpart 
and one-fourth have radio detection. 

We note here that half of our sample overlaps the spec- 
troscopic sample discussed by Treu et al. (2005). 

The present paper is structured as follows. Section 2 
describes the multi-wavelength data-set and Section 3 our 
galaxy sample selection. In Section 4 we give details about 
our data analysis, by exploiting in particular the SEDs of 
our galaxies, the information from multi-colours and broad- 
band flux ratios. Section 5 discusses and critically analyses 
these results, while Section 6 summarizes our conclusions. 

We adopt a cosmological model with Q m —0.3, Oa=0.7, 
and # = 70 km s" 1 Mpc' 1 . 



2 THE DATA 

With the aim of studying the multiwavelength photometry 
of extragalactic sources measured by different instruments, 
we need to measure the bulk of the emission from each object 
in each photometric band. Only with such kind of approach 
an SED can be considered reliable. In our work, we per- 
formed aperture photometry in most bands, and decided to 
bring the measurements as close to the total magnitudes as 
possible (see also Yan et al. 2005). 



2.1 Optical Imaging 

The GOODS project (Dickinson et al. 2003) has sur- 
veyed with ACS/HST two separate fields, the Chandra 
Deep Field South (CDFS) and the Hubble Deep Field 
North (HDFN) with four broad-band filters: F435W (B), 
F606W(V), F775W(i), and F850LP(z). The analysis pre- 
sented here is based on the version 1.0 of the reduced, 
stacked and mosaiced images for all the data acquired over 
the five epochs of observation (Giavalisco et al. 2004). 

Source extraction and photometric measurements in 
the B, V, i and z bands have been performed by run- 
ning a modified version of SExtractor (Bertin & Arnouts, 
1996). We have considered the total magnitudes (as from 
the BEST_MAG output parameter in SExtractor, see Gi- 
avalisco et al. 2004). 

Moreover, optical and near-IR ground-based photomet- 
ric imaging in several bands for the GOODS-North field was 
taken from Capak et al. (2004). The [/-band data were col- 
lected using the Kitt Peak National Observatory (KPNO) 
4m telescope with the MOSAIC prime focus camera. The B-, 
V-, R-, I-, and z'-band data were collected using the Sub- 
aru 8.2 m telescope and Suprime-Cam instrument (Miyazaki 
et al. 2002). The HK' data were collected using the QUIRC 
camera on the University of Hawaii 2.2 m telescope (Hodapp 
et al. 1996). Details on the data reduction and photometric 
analysis are reported in Capak et al. (2004). A 3" diameter 
aperture was used for the photometry. We applied the cor- 
responding aperture corrections reported by Capak et al. in 
all bands. 



2.2 Morphological Analysis 

Our basic selection is a sample of morphologically classified 
elliptical galaxies in the GOODS-N based on the imaging 
data. We used the publicly available catalogue of Bundy et 
al. (20050, 

which was constructed with a magnitude limit of 
zab = 22.5 over an area of 160 arcmin 2 . The morphological 
selection is based on the visual inspection of ACS version 
1.0 data. Our sample selection is based on a catalogue of 
168 bona-Sde E/SOs (classes and 1 in Bundy's catalogue) 
out of a total sample of 1576 galaxies brighter than this 
magnitude limit. 

2.3 Deep Spitzer Imaging in the Near- and Mid-IR 

As part of the GOODS project, the Spitzer Space Telescope 
has surveyed the GOODS-N field in the IR between 3.6 and 
8.0 ^m using IRAC and in the range 24-160 ^m with MIPS. 
In this paper we have made use of the public reduced data 
released by the GOODS team and available on the Wet0. 

For the IRAC data, we have adopted the official 
GOODS reduced images (epoch 1 + 2) and performed our 
own photometry. The IRAC source identification has been 
performed with SExtractor independently in the four maps 
at 3.6, 4.5, 5.8 and 8.0 urn. 

In order to obtain the most accurate Spectral Energy 
Distributions (SEDs), we have computed the fluxes of each 
source by performing aperture photometry in the four IRAC 
bands at the positions originally detected in the ACS z- 
band. Assuming that essentially all the sample sources are 
seen as point-like by the IRAC ~ 2 arcsec FWHM PSF 
imager, we computed with SExtractor (Bertin & Arnouts, 
1996) the fluxes within a 3.8 arcsec diameter aperture. This 
choice is supported by an accurate analysis performed by the 
SWIRE teairfl They constructed color-magnitude diagrams 
for various types of objects, in particular main-sequence 
stars. It was then found that the scatter in these diagrams 
is minimized through the use of an aperture of 3.8 diam- 
eter aperture, and corresponds to roughly twice the beam- 
width. To obtain total fluxes, we then applied the correction 
factors indicated by the SWIRE tearrjj. We have indepen- 
dently verified that the IRAC/SWIRE aperture corrections 
are consistent with those derived by fitting the radial bright- 
ness profiles of few stars in the GOODS fields (bright stars 
where selected on the basis of their visual morphology in 
the ACS bands). For the only one source spatially resolved 
(# 15), we used Kron like magnitudes (AUTOJV1AG output 
parameter in SExtractor). For three blended IRAC sources 
present in our sample, we applied a deconvolution procedure, 
as described in Section 3.1. 

The public MIPS data-set includes calibrated maps and 
a catalogue of 24 um sources with flux densities S_4 > 80^iJy. 
The public photometry is based on a PSF fitting algorithm, 
where the SExtractor positions of the IRAC sources are used 
as input to the MIPS source extraction process. The MIPS 
24 am PSF was generated from isolated sources in the im- 
age, and re-normalized based on the aperture corrections 

1 available at http://www.astro.caltech.edu/GOODS_morphs 

2 http:/ /data. spitzer. caltcch.edu/popular/goods 

3 http: / / data.spitzer.caltech.edu/popular /swire/20050603_enhanced_ 

4 See Note 2 
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Figure 1. Top panel: optical images of the 16 spheroids with MIR emission and reliable optical associations. Bottom panel: optical 
images of the 14 spheroids with X-ray emission and undetected in the MIR. For each source a cutout of the z' band ACS images is 
presented. The identification numbers are those reported in Table 3. Each postage has a size of 8"x8". North is up, East is at left. 



published in the MIPS Data Handbook (v2.1, section 3.7.5, 
table 3.12). 

IRS Peak-Up imaging data at 16 fj,m are also available 
from Teplitz et al. (2005) over an area of 35 square arcmin- 
utes in the GOODS-N to an average 3 sigma depth of 0.075 
mjy, for a total 149 detected sources. Three of these are 
in common with our spheroidal galaxy sample. The 16 /im 
photometric band provides unique information on the shape 
of the observed SEDs. In fact, the combination of the com- 
plete Spitzer filters set represents a powerful tool to sample 
the evolution of the PAH features. For example, the 7.7 fim 
PAH feature can be traced from the local universe, where it 
lies in the IRAC channel 4 passband, to 0.8 < z < 1.3, where 
it falls in the 16 fim IRS Peak-Up filter, up to redshifts near 



z < 2, where it enters the 24 jj,m MIPS filter (Teplitz et al. 
2005). 



2.4 X-ray data 

Very deep X-ray observations (2 Megaseconds in total) of 
the Chandra Deep Field North (CDFN) have been reported 
by Alexander et al. (2003) . We used them to search for nu- 
clear activity in spheroidal galaxies of the GOODS-North 
field in addition to those obtained from the IR emissions. 
The X-ray data have been retrieved from the public archive 
and processed with standard tools making use of the cali- 
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Figure 2. Overlays of the MIPS 24 (im contours over the op- 
tical ACS i' band images for spheroidal galaxies with multiple 
IDs inside the 24/im error-box (blended candidates). From up- 
per left panel to right panel ID #0 (z=0.41), ID 4 (z=0.845), ID 
11 (z=0.766). The lower panel corresponds to overlays with the 
IRAC 3.6/mi contours. Each postage has a size of 8"x8". North 
is up, East is to the left. 



brations associated with the CIAC0 software (version 3.2.1). 
Each of the 20 pointings of the Chandra dataset has been 
registered on a reference pointing (OBSID 3293) and aligned 
using the align_ev 

10 

tool. The charge transfer inefficiency 
and gain corrections were also applied to each single point- 
ing. X-ray images were accumulated in the soft (0.5-2.0 
keV), hard (2-8 keV) and full (0.5-8 keV) bands. 

X-ray counts in the three bands were extracted within 
a circular region whose radius is a function of the off-axis 
angle (typically of the order of 2-3 arcsec). The background 
was estimated locally for each source as the average of counts 
in several source-free surrounding regions. The number of 
counts and the X-ray fluxes are fully consistent with the 
values published in Alexander et al. (2003) , and are reported 
in Table 3. 



2.5 Optical spectroscopy 

In the last few years various observational programs have 
undertaken a systematic spectroscopic follow-up in the 
GOODS-N. In particular, the Keck team (Wirth et al. 2004) 
reported the results of an extensive imaging and spectro- 
scopic survey in the GOODS-North field with DEIMOS on 
the Keck-II telescope. Observations of 2018 targets from a 
magnitude-limited sample of 2911 objects to Rab ~ 24.4 
yielded secure redshifts for a sample of 1440 galaxies and 
AGNs, plus 96 stars. All the Keck spectra are publicly avail- 
ably All our objects have a redshift from the Wirth et al. 
(2004) catalogue. However, we did not found the correspond- 
ing spectrum in the Keck archive for 10 sources in our sample 
(IDs 3, 4, 9, 10, 15, 16 and IDs 3x, 4x, 8x and lOx as in Table 
3). We then recovered four of the missing spectra (IDs 4, 9, 
10 and 16) from the sample of Treu et al. (2005). All sources 



have a z-quality value of 4 (99% confidence level) or 3 (90% 
confidence level), as in the code assigned by Keck team. 



3 SAMPLE SELECTION 

3.1 Spheroids with mid-IR excess emission 

As a first step, we cross-correlated the optical catalogue in- 
cluding the 168 morphologically selected spheroids (see Sec- 
tion [272]) , with the bright Spitzer — MIPS 24 /jm catalogue 
provided by the GOODS team (S24 > 80^Jy, Section [23)) . 
For 21 of the 168 spheroids we found a MIR counterpart 
within a distance of 2" from the optical position. We checked 
the morphology of each of the 21 candidates taken from the 
Bundy et al. (2005) catalogue on the i-band ACS image. We 
then excluded two objects (#14 and #18) from the sample 
because they look like late- type galaxies rather than E/S0 
to our visual inspection (however, we have verified that the 
inclusion of these two sources does not change the main con- 
clusions of this work). 

Source confusion is a potential problem for the optical 
identification of the 19 MIR sources in the final sample: the 
very different angular resolution of the Spitzer and HST 
instruments enhances the probability of finding more than 
one ACS source falling within a MIPS beam. However, for 
the majority of our sources the optical to NIR and MIR 
identifications turned out to be straightforward. Only six 
of the optically-selected galaxies present extended IRAC or 
MIPS emissions and confusion problems. For three of them 
(#9, #10 and #16) we could use the spatial information 
available in the HK'-band image to deconvolve the 3.6-8.0 
(im fluxes (for details on this procedure see Franceschini et 
al. 2006). Using the HK' positions for reference, we applied 
a PSF-fitting algorithm based on IDL procedure written ad- 
hoc to deblend also the 24 fj,m emission. 

The confusion issue for the three remaining objects (ID 
#0, #4 and #11) cannot be resolved by just looking at the 
HK' image. We will consider these 3 blended sources sepa- 
rately, and their MIR fluxes as upper limits of the spheroidal 
galaxy emission, while the final catalogue of robust unam- 
biguous associations includes 16 E/SOs with 24 /jm emission 
above S24 > 80^iJy. 

Figure [1] (upper panel) shows optical postage stamps of 
the 16 reliable spheroids in our final sample. For each source 
a cutout of the four ACS images (B, V, i' and z' bands) is 
presented. The identification numbers are those appearing in 
Table 3 (see also £|3,3|l . Postage stamps of the three blended 
sources are reported in Figure _I1 where the optical ACS i' 
band images are overlaid on the MIPS 24 /im contours. For 
the 16 isolated spheroids, the association of the 24 fim and 
the optical/near-IR source is unambiguous. 

For what concerns the sources affected by confusion 
problems, the X-ray counterparts, where available, (see Fig- 
ure below) support the hypothesis that the spheroidal ob- 
ject is strictly correlated with the stronger IR emitter, under 
the assumption that the two emissions are related. 



5 http://cxc.harvard.edu/ciao/ 

6 http: //cxc. harvard.edu/cal/ASPECT/align_evt/ 

7 http:/ /www2. keck. hawaii.edu/science/tksurvey 



3.2 Spheroids with X-ray emission 

More than half (12/19) of our Spitzer-MlPS sources have a 
unique X-ray counterpart within a search radius of 2 arcsec. 
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Figure 3. Overlays of the i-band images with X-ray Chandra contours for 4 spheroidal galaxies with confused identifications for the 
MIPS 24 source. 



Four of the 12 X-ray sources are associated with the objects 
showing confusion problems at 24 /im (see § 13.10 . By exploit- 
ing the arcsec X-ray image quality, the soft X-ray contours 
of the 4 Chandra sources were overlaid on the HST/ACS 
images (Figure [3J. The X-ray/HST association is obvious 
for 3 of them where the X-ray emission is centered on the 
spheroidal component. The closest counterpart of source 10 
is also a spheroidal galaxy, although the correct association 
is not so clear. 

For six out of the seven MIR spheroidal galaxies un- 
detected in both soft and hard X-ray bands, upper limits 
were computed by X-ray counts in circular regions (2 arc- 
sec of radius) of the CDFN map centered and reported in 
Table 3. The optical position of source 19 was too close to 
a bright X-ray source to obtain a reliable estimate of the 
X-ray upper limit. 

In addition, we found 14 associations of Chandra X-ray 
sources among the 147 z-band selected spheroidal galaxies 
in the GOODS-North area without detectable MIR emis- 
sion. The optical z-band images for all of them are reported 
in Figure [T] (lower panel). All appear to display standard 
elliptical morphologies. Photometric and other data for this 
sub-sample are reported in the third panel of Table 3. The 
majority of these sources (8 out of 14) are detected in hard 
band (2-8 keV). 

To summarize, these 14 additional X-ray luminous ob- 
jects bring our total sample of spheroidal galaxies in the 160 
sq. arcmin GOODS-North field with signs of activity (either 



from MIR or X-ray data) to a total of 33 out of 168 objects 
in the original optical sample (corresponding to a fraction 
of 20%). 



3.3 Summary of the multiwavelength photometry 
for spheroids in the GOODS-N 

We have built a multi- wavelength photometric catalogue for 
our sample of spheroidal galaxies with MIR or X-ray emis- 
sions in the following way. Starting from the optical z-band 
positions, we used a search radius of 1" in each bands to look 
for the cross-identifications in the corresponding catalogues. 
As mentioned in Section 13.11 direct visual inspection of the 
ACS images guarantees the reliability of the final counter- 
parts. 

The photometry for each source is presented in Table 
3. We adopt a common value of 10% (15%) of the measured 
fluxes as photometric errors for the IRAC (MIPS) bands, 
in order to reflect the systematic uncertainties of the instru- 
ments The main contributions to these uncertainties are due 
to the colour-dependence in the flat field and to the absolute 
calibration (see for example Lacy et al., 2005, and the IRAC 
and MIPS Data Handbook). 

For what concerns the optical photometry, in the SED 
fitting analysis we adopted for each source the U, B, V, R, I 
and z' band from the Capak et al. (2004) paper (Section l2.ll) . 
We used the ACS photometry only when this information is 
missing (see Table 3). 
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All the 19 sources in the final catalogue of spheroids 
with 5*24 > 80/iJy and the 14 additional X-ray detections 
only are found within the spectroscopic Keck catalogue (Sec- 
tion 12. 5[) . and spectroscopic redshifts are also reported in 
Table 3. Table 3 is splitted in three sections, one referring 
to the catalogue of 16 MIPS 24/im sources with unambigu- 
ous optical identifications, the second to the 3 blended MIR 
sources, and the last one to the 14 spheroidal galaxies with 
X-ray detection only. For the potentially confused sources, 
Spitzer (IRAC+MIPS) fluxes have to be considered as up- 
per limits only. 



4 DATA ANALYSIS 

In this Section we exploit the previously discussed multi- 
wavelength data to test the origin of the observed emissions 
and verify if these may be ascribed to the standard pas- 
sively evolving stellar populations, or instead require star- 
formation or AGN activities to take place in them. 

To this end, we use the optical-IR Spectral Energy Dis- 
tributions (SEDs) to check for the presence of hot dust emis- 
sions by a nuclear power-source (AGN) or emission by warm 
dust in star-forming regions. We will also look for the pres- 
ence of emission lines the optical spectra and check for the 
ionization level and line width to discriminate AGN activity. 

Then optical to IR colours, X-ray to optical and X-ray 
to IR flux ratios are used as diagnostics of the emission pro- 
cesses. In a few cases the very deep X-ray data will provide 
unique diagnostic information. 



4.1 Analysis of the Spectral Energy Distributions 
for the IR-emitting Spheroids 

In order to understand the nature of dust emission by 
the IR-emitting objects, we have assembled in Figures 
H and [H the SEDs from UV to MIR for all of them. 
We have first attempted to interpret these observational 
SEDs of our spheroids galaxies using an automated fit- 
ting routine contained in the "Le Phare" code (available 
at http : //www. lam. oamp. fr/arnout s/LE_PHARE. html ). 

Le PHARE (PHotometric Analysis for Redshift Estima- 
tions, Ilbert et al. 2006) is a publicly available set of FOR- 
TRAN programs aiming at computing photometric redshifts 
through the best-fitting SED analysis. Le PHARE's struc- 
ture is based on two parts: 

• A preparation phase, composed of different programs 
in order to select theoretical SEDs, filters and build theoreti- 
cal magnitudes. The different programs have been developed 
in order to easily get basic information relative to the filters 
(A me an, AB-corrections ...) and the SEDs (k-correction ver- 
sus z, color-color diagrams for stars, quasars and galaxies 
...). 

• A running phase, based on the photometric redshift 
program itself. The program is based on a simple x 2 fitting 
method between the theoretical magnitudes and an observed 
photometric catalogue. 

For our galaxies we already know the spectroscopic red- 
shift, so we have fixed z and have used Le PHARE only 
to obtain the best-fitting SED through a comparison with 
21 templates of galaxies (including Ellipticals of different 



ages, lenticulars, Spirals and Starbursts) and AGNs (includ- 
ing low- and high-ionization QSOs, reddened QSO, Seyfert 
1, 1.8 and 2 and two kinds of ULIRGs) provided by M. Pol- 
letta et al. (in preparation). 

Since the SED-fitting is based on the x 2 method, the 
most relevant photometric points are obviously those with 
smaller photometric errors (i.e. typically the optical/NIR 
ones) . 

From this first exploratory SED fitting we have inferred 
a rough spectral classification, and have looked for evidences 
of a MIR excess in the observed 24 /im fluxes compared with 
the best-fitting model spectra. 

The left panels in Figures [3] and [S] provide a comparison 
of our model fits based on local templates with the obser- 
vational SEDs for the spheroidal sample galaxies. For each 
source we report in the plot a reference indicating the na- 
ture of the best-fit model derived by the Le Phare code and 
the reduced x 2 corresponding to that solution. The large 
values of the \ 2 statistics are mainly caused by the inability 
of the adopted SED templates to reproduce the observed 24 
fim emission (but also the 5-8 fim fluxes of various sources) . 
We note that the primary solution of this photometric anal- 
ysis is quite robust, within the uncertainties. In fact, only 
for a couple of sources the secondary solutions derived by 
Le Phare are acceptable (in terms of minimization of the 
X 2 ) and imply a different spectral type (from an SED domi- 
nated by star formation in favour of and SED dominated by 
an AGN, or viceversally) . However, the main implication of 
this analysis is that none of our sample spheroidal objects 
is reproduced by a typical old and passive stellar popula- 
tion. On the contrary, only late-type star-forming spirals or 
Seyfert galaxies seem to better match the observed SEDs. 
This already implies that our selection naturally favours the 
detection of a peculiar class of active ellipticals (see discus- 
sion below). 

We see that in many cases (6 out of 19) the observed 
24 fim fluxes show a significant excess compared with the 
expectations based on the local spectral templates matched 
to the optical- NIR data. We discuss here possible interpreta- 
tions of such MIR emission. Three different hypotheses are 
considered: 1) mass loss from late-type stars (Knapp et al. 
1992), 2) IR emission by warm dust and PAH molecules re- 
lated with star-forming regions (Madden et al. 1999; Xilouris 
et al. 2004), 3) dust emission by circum-nuclear torii heated 
by an AGN (Knapp et al. 1999). 

The first interpretation follows from the study by 
Knapp et al. (1992) of a sample of nearby elliptical galaxies 
detected in the MIR (12 fim) by the IRAS surveys. The au- 
thors, after considering several possible origins for the 12 fim 
emission, conclude that for almost all the normal non-AGN 
ellipticals the MIR emission is likely due to photospheric 
and circumstellar emission by evolved red giants stars in 
the galaxies. The MIR (12 fim) and NIR (2.2 (im) fluxes of 
the nearby ellipticals are compared by Knapp et al. (1992) 
to those of evolved red giants, showing the existence of a 
large population of mass loosing stars. MIR and NIR fluxes 
are then used by them to derive rates of stellar mass loss. 

To test the hypothesis that the MIR emission from 
early-type galaxies in our sample has the same origin of the 
nearby ellipticals, we have compared the NIR and MIR prop- 
erties of the two samples. Since our galaxies are at a median 
redshift < z >~ 0.68, from our observed 3.6 and 24 /im 
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Figure 4. Spectral Energy Distributions of spheroidal galaxies with detected 24 /im flux. The left panels show a comparison of local 
template spectra with the observational SEDs. The right panels illustrate two-component fits: an evolved stellar population (dashed line) 
and dust obscured AGN emission (green line). The sum of the two components is shown as the solid red line. 
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Figure [4] (continued) 
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Figure 5. Same as for Figure[4]for the three MIR sources with multiple optical associations (blended sources, see images in Figure[2]|. 
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fluxes we have computed rest-frame luminosities at 2.2 and 
12 /im with suitable K-corrections (that turned out to be 
rather insensitive to the adopted template spectrum), and 
then compared them in Figure [6] with intrinsic 12 and 2.2 
/im luminosities of local objects. 

The values of the MIR to NIR luminosity ratios for our 
early-type galaxies (open stars) turn out to be about a factor 
of ten higher on average than for the Knapp et al. nearby el- 
lipticals (filled circles). This would imply a fraction of 12 /jm 
luminosity produced by mass loss (and the relative mass loss 
rate) in our galaxies far too high if compared with the values 
found for local galaxies with MIR excess. We conclude that, 
if some contribution to the observed MIR emission by mass 
loss from evolved giant stars is present in our spheroidal 
galaxies, it is certainly not the dominant energy source at 
these wavelengths. 

The presence of diffuse dust emission, responsible for 
the observed MIR excess in the spectra of early-type galax- 
ies, has been considered by several authors (i.e. Madden et 
al. 1999; Xilouris et al. 2004). To test this hypothesis for 
our galaxies, we have considered the results of the automatic 
SED fitting described in the previous section and shown in 
the left panels of Figs. [4] and [5] In 8 galaxies out of 19, 
the optical-to-MIR SEDs can easily be fitted with normal 
galaxy templates and residual star-formation. In at least 6 
further cases there is instead evidence that the rest-frame 
MIR emission is significantly in excess of any reasonable 
level that might be explained by diffuse dust in star-forming 
regions. 

Our third considered hypothesis is therefore that the 
MIR excess observed in the GOODS-N early-type galaxies 
is contributed by the presence of AGN circum-nuclear dust 
re-radiation. This interpretation is suggested by the shapes 
of the observed IR SEDs, which exclude the alternative pos- 
sibility of direct non-thermal quasar unobscured emission. 

We have then tried to fit the observational data with 
a combination of two different spectral components: an 
evolved stellar population to reproduce the optical and NIR 
spectrum, plus an obscured AGN emission component to 
explain the MIR data. The synthetic SEDs of the old com- 
ponent, including a mix of passively evolving stellar pop- 
ulations and younger stars, have been modelled with a 
population-synthesis code reported by Berta et al. (2004), 
based on the Padova stellar isochrones. The synthetic spec- 
tra are taken from a set of templates used in Fritz et al. 
(2006, submitted). We found that sources ID5 and ID6 are 
reproduced with a typical evolved old stellar population of 
1 Gyrs, while ID 10 and ID 5 with an older population of 4 
Gyrs. Other six sources (IDs 0, 1, 3, 8, 9 and 17) are better 
represented by a galaxy during a post-starburst phase with 
a second main episode of star formation -forming ~ 10% of 
the total stellar mass- in the range ~ 10 8 — 3 x 10 9 years. 
The remaining objects (IDs 2, 4, 7, 12, 13, 15, 16, 19 and 20) 
are instead fitted with a main-impulsive burst and a con- 
tinuous star formation which is truncated at 2 x 10 7 years. 
The best-fit solutions for the stellar component are also used 
to provide a measure of the barionic mass of each source (a 
Salpeter IMF in the range 0.15-120 Mq has been consid- 
ered). The stellar masses are reported in the last column of 
Table 1, with a median value of log(M) ~ 10.7 Mq. 

The AGN IR component is modelled as the emission by 
a dusty torus heated by a central AGN, as detailed below. 
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Figure 6. Relation of the 12 /im versus 2.2 /im rest-frame lu- 
minosities for MIR emitting spheroidal galaxies. The GOODS-N 
early-type galaxy colours (open circles) are compared to those 
of nearby elliptical galaxies studied by Knapp et al. (1992; filled 
circles), and show a large excess in the 12 fim luminosity. 



The right-hand panels in Figs. [4] and [5] illustrate the two 
component (old population + dusty torus) solutions. 

Our adopted torus model has been discussed in Fritz et 
al. (2005), to which we refer for a more detailed explanation. 

The model reproduces the AGN spectrum as a combi- 
nation of a nuclear non-thermal spectrum plus a dusty torus 
emission. A flared disc geometry, i.e. two concentric spheres 
with the polar cones removed, is adopted for the latter. The 
torus size is defined by the angular opening angle and by the 
radius of the outer sphere, while the inner radius is set by 
the sublimation temperature of graphite and silicate grains 
under the influence of the strong nuclear radiation field. 

The computation of the thermal and scattering emis- 
sions throughout the torus is performed by numerically solv- 
ing the radiative transfer equation with an iterative proce- 
dure. 

As discussed in Fritz et al., the model provides excellent 
fits to the observed SEDs of a wide variety of active nuclei 
and quasars. 



In almost all cases, the combination of an evolved stel- 
lar population fitting the optical and NIR spectrum, and 
an obscured AGN emission component to explain the MIR 
spectrum, always provides much better fits to the data (as 
can be argued from the comparison of the \ 2 statistics in the 
left and right panel of Figures [4] and dj. As anticipated, in 6 
cases this solution including an AGN is required by the IR 
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Figure 7. Observed optical spectra of spheroidal galaxies with excess MIR emission. The ID numbers are those reported in the top 
panel of Table 3. The dashed region indicates an absorption feature of the sky. 



12 G. Rodighiero et at 




wavelength (angstroms) 



Figure [7] (continued) 
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Table 1. Object Classification 



ID 


X-ID 


z 


optical spectral 


Overall SED 


Mass 








classification 


classification 


(M ) 


1) 




0.410 


ELG 


AGN 


10.88 


1 


X-067 


0.038 




SB 


12.15 


2 


X-082 


0.679 


AGN2 


AGN 


11.21 


3 


X-110 


1.141 




AGN 


11.34 


4 


X-113 


0.845 


LINER 


AGN 


11.32 


5 




0.156 


ELG 


SB 


10.05 


6 




0.512 


LINER 


AGN? 


9.50 


7 


X-115 


0.680 


AGN2 


AGN 


10.65 


8 


X-149 


1.223 


ELG 


SB/AGN? 


10.78 


9 


X-160 


0.848 


Starforming 


AGN 


11.50 


10 


X-169 


0.845 


Early- type 


AGN? 


11.10 


11 




0.766 


EarIy+[OII] 


SB 


11.02 


12 




0.277 


Early- type 


SB 


9.76 


13 


X-240 


0.961 


AGN2 


AGN 


12.50 


15 


X-383 


0.105 




SB 


10.25 


16 


X-388 


0.231 


Early-type? 


SB 


10.40 


IT 




0.899 


Starforming 


SB 


10.20 


19 




0.278 


Starforming 


SB 


9.44 


20 




0.911 


Early+faint[OII] 


AGN 


10.74 


lx 


X-414 


0.798 


Early 


SB 


10.90 


2x 


X-230 


1.013 


Early-type, post-SB 


SB? 


11.51 


3x 


X-294 


0.474 




AGN 


10.14 


4x 


X-249 


0.475 




SB 


10.56 


5x 


X-286 


0.954 


Early 


AGN 


11.40 


6x 


X-257 


0.089 


Early- type 


Inactive 


9.79 


7x 


X-274 


0.321 


Early-type 


Faint AGN? 


10.74 


8x 


X-354 


0.569 




SB? 


10.44 


9x 


X-241 


0.850 


Early-type 


AGN 


10.66 


lOx 


X-210 


0.848 




SB? 


11.26 


llx 


X-212 


0.943 


Early-type 


AGN 


10.90 


12x 


X-194 


0.556 


Early+faint[OII] 


AGN 


10.34 


13x 


X-131 


0.632 


Starforming 


SB 


10.20 


14x 


X-114 


0.534 


Early-type 


SB 


10.66 



data. (Obviously, we did not attempt to constrain the torus 
physical parameters with the few available data-points). 

In order to have a complete census of the physical prop- 
erties of the spheroidal population, we checked the observed 
optical-to-IR SEDs of the z-band selected galaxies lacking 
a detectable mid-IR or X-ray counterpart (135 out of the 
initial 168). For these sources we built a multiwavelength 
photometric data-set similar to that described for the main 
sample of this work, covering the spectral range 0.3-8.0 /im. 
The vast majority of these objects show SEDs consistent 
with that of a passively evolving stellar population, typi- 
cal of quiescent elliptical galaxies. We found that only few 
sources (~6%) present some indication of "activity" in the 
observed SEDs, simply because the photometric data points 
at wavelengths larger than 5 /xm show an underlying rising 
spectra. However, a more detailed work on this photometric 
data-set would be required to constraint the properties of 
the full morphological spheroidal sample. 



4.2 Optical Spectroscopy 

The optical spectra available for our high-redshift spheroidal 
galaxies are reported in Figures [7] and [8] The emission lines 
were analysed only in terms of the equivalent widths and 
flux ratios of contiguous lines, because the optical spectra 
retrieved from the WEB were not flux calibrated. Objects 
with a red continuum and faint (or absent) emission lines 
(EW < 5A) are classified as passive early-type galaxies. 
When two or more emission lines were measured, we applied 
the usual diagnostic diagrams (Veilleux & Osterbrock, 1987) 



in order to classify the galaxies according to the star-forming 
emission line, Seyfert-2 or LINER categories. In few cases we 
could not disentangle the Seyfert-2s/LINERs classes due to 
the limited observed spectral range. Sources with more than 
one emission line were classified as Starforming, while when 
only a single emission line is visible, we classify the object 
as an Emission Line Galaxy (ELG). 

Optical spectra for our sub-sample of galaxies detected 
only in the X-ray bands are reported in FigureJU and have 
been analysed based on similar procedures. 

Since the original parent sample is morphologically se- 
lected spheroidal galaxies, we might be surprised by the high 
incidence of strong emission lines in the optical spectra (67% 
of the objects show at least one line with EW > 9 A), con- 
sidering the usual fairly good agreement between the mor- 
phological and spectroscopic classifications. However, this 
fraction of optical emission line spectra is not untypical if 
referred to the original parent sample: a 10% of emission- 
line objects among high-z morphological ellipticals was also 
found in the K20/GOODS sample by Cassata et al. (2005). 
Moreover, it is well known that the fraction of spheroids with 
emission lines increases with redshift, with a typical value 
of ~20% in the redshift range 0.5 < z < 1 (Brinchmann et 
al. 1998; Treu et al. 2002; Treu et al. 2005). In particular, 
~25% of the spheroidal galaxies in the Treu et al. sample, 
when limited to the same morphological classes considered 
in this work (class=0,l), present Oil emission lines (~10% 
with EW > 5 A). This might suggest that the fraction of 
active spheroids based only on optical indicators is consis- 
tent with that determined with the combination of X-ray 
and MIR selections. 



4.3 X-Ray Properties 

We consider here the X-ray emission properties of spheroidal 
galaxies. A first indication comes from Figure [5] plotting the 
soft and hard X-ray fluxes versus the R-band magnitudes. 
The first point to be remarked is that the optical and X-ray 
fluxes appear largely unrelated. 

The values of the X-ray-to-optical flux ratio (X/O = 



log if = log 



fx + C + with C=5.5) can then be used 



to obtain a preliminary classification of the X-ray emission 
components (e.g. Maccaccaro et al. 1988): various classes of 
X-ray emitters are characterized by different values of X- 
ray-to-optical flux ratio. Optically and X-ray selected AGN 
have, on average, relatively high X/O values (-1< X/O <1), 
while infrared emitting starforming galaxies show lower X/O 
ratios (-2 < X/O < -1; Alexander et al. 2001; Frances- 
chini et al. 2003; Ranalli et al. 2003), and even lower values 
(X/O < -2) are typical for X-ray quiescent normal galaxies, 
whose X-ray flux is due to the integrated emission of com- 
pact sources and /or hot interstellar plasma (Hornschemeier 
et al. 2001). 

This classification based on the X-ray-to-optical flux 
ratio should be considered with caution and is somewhat 
dependent on the adopted energy range. For this reason we 
consider separately the 0.5-2 keV and 2-8 keV X-ray fluxes 
in the figure. The six sources detected in the hard band 
(filled squares in bottom Figure [9]) are well within the AGN 
locus. Their soft X-ray fluxes are lower than the 2-8 keV 
ones and the observed 0.5-10 keV luminosities are within 
the Seyfert range (10 42-44 erg/s). 
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Figure 8. Optical spectra of spheroidal galaxies with excess X-ray but no MIR emission. The ID numbers are those reported in the 
bottom panel of Table 3. 



Three of them (IDs #2, #9, #20) fall outside the AGN 
region ( — 1 < X/O < 1) in the soft band diagram, and, ac- 
cording to the classification based on the X-ray-to-optical 
flux ratio, they would be classified as star-forming galax- 
ies. However, the AGN locus boundaries in X/O diagrams 
should not be considered as strict limits. From an inspection 
of the two panels of Figure 9, it is clear that the three sources 
have hard spectrum' This hypothesis is corroborated by the 
results of the spectral analysis. 

Five out of the six hard X-ray detected sources (2, 4, 
7, 9, 13) have enough counts (> 100) to allow us to per- 



form a moderate-quality spectral analysis. Source and back- 
ground spectra, response matrices and effective areas have 
been extracted from the merged events file using the stan- 
dard CIAO tools developed to properly weight responses and 
effective area files for multiple extraction regions (mkwarf 
and mkwrmf ; see Civano et al. 2005 for a detailed descrip- 
tion). Source spectra were re-binned in order to have enough 
counts per bin (15-20) to apply \ 2 statistic, and were fitted 
with XSPEC 11.3 (Arnaud 1996). Power-law model spectra 
with intrinsic absorption at the source redshift have been 
adopted. The results of spectral fitting are reported in Ta- 
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Figure 9. Top panel: the R-band magnitude versus the 0.5-2 keV X-ray flux for the 12 Spitzer-MlPS ellipticals with X-ray counterparts 
(filled squares), those without X-ray counterparts (open squares), and the spheroids detected only in the X-rays but not in the MIR 
(crosses). The area between the two continuous lines represents the locus typically occupied by known AGNs, and corresponding to 
log(XfO) = ± 1. Bottom panel: the same plot against the 2—8 keV X-ray band. 



ble 2. Errors are at the 90 per cent confidence level for one 
interesting parameter (A^ 2 = 2.71). The best fit Nh values 
and the corresponding unabsorbed 0.5-8 keV luminosities 
(10 43 -10 44 ) are typical of moderately obscured Compton- 
thin Seyfert galaxies (Risaliti et al. 2002). 

Despite the relatively low counting statistics, the X-ray 
spectrum of ID #9 (Figure fTTJ|) does clearly require a rather 
large column of cold gas {log{Nu) > 23 cm -2 ). Source #2 
has a very flat spectrum and an emission line at the redshift 
expected from the iron K a with a large (1 keV) EW is ev- 
ident from the residuals of a power law fit (Figure 1 10 p . A 
similar model provides a good description of the ID #20's 
spectrum, though with lower significance due to the poor 
counting statistics. In both cases, Compton thick absorp- 
tion is the most likely explanation of the X-ray properties. 
Assuming that the observed luminosity is entirely due to 
reflection by thick matter with an albedo of 2% (Comas- 
tri 2004), the intrinsic luminosity would be of the order of 



Table 2. X-ray spectral analysis parameters 



AID 


id 


z 


r 


NHxlO 22 


chi 2 /dof 


L x (0.5-10 k 




2 


0.679 


0.40±0.33 






O.lxlO 43 


113 


4 


0.84 


1 56+ ' 12 


<0.6 


39/31 


l.OxlO 43 








1.8 


0.7±0.3 


41/32 


l.OxlO 43 


11B 


7 


0.68 


1.57±0.07 


49+ ' 22 
u -0.13 


63/62 


2.1X10 43 








1.8 


0.84±0.13 


73/63 


2.1X10 43 


160 


9 


0.848 


i.iot°;|° 


— 5.05 


6/11 


8.7X10 42 








1.8 




8/12 


9.7X10 42 


240 


13 


0.961 


1.69±0.09 




105/96 


1.3X10 44 








1.8 


5.59±0.36 


112/97 


1.6 X10 44 



a Source number 82 has been fitted with the Cash statistic. 

about 10 44 erg/s for source #2 and about 10 43 erg/s for 
source #20. 

Finally source #3 is also classified as an AGN on the 
basis of the observed X/O in the soft band (it is undetected 
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Figure 10. The X-ray spectra of sources 2,9 and 20 (from top to 
bottom)/ An absorbed power law model with both spectral index 
and absorption free to vary has been adopted for source #9. The 
flat X— ray spectrum of source #2 is well fitted by a single power 
law model (see Table 2). The low counting statistic of source #20 
does not allow to perform a proper spectral fit. A power law with 
slope fixed at T = 1.8 is adopted to show excess emission around 
the energy expected for a redshifted iron Ko line. 



in the hard band) and the relatively high (~ 7 x 10 42 erg/s) 
X-ray luminosity. Therefore, at least 7 of the 12 X-ray de- 
tected sources turn out to be luminous AGN. 



4.4 The sample of X-ray detected active 

spheroidal galaxies with no MIR emission 

We have attempted to complement the information provided 
by the limited sensitivity of the 24 /im map by including in 
our analysis also the sample of 14 spheroidal galaxies de- 
tected in X-rays but not in the MIR. In these cases, ongoing 
activity is traced by the excess X-ray emission. Detailed in- 
formation on these objects is reported in the third panel of 
Table 3. 

The crosses in Figure[9]indicate the soft and hard X-ray 
fluxes against R magnitudes for such X-ray detected objects 
(numbers refer to the ID's reported in the third panel of Ta- 
ble 3). The sources occupy a region of the plot intermediate 
between that of AGNs and that of starbursts. For 5 objects, 
those with ID #3x, 5x, 9x, llx, 12x, an AGN predominance 
in the X-ray flux is suggested. 

The X-ray luminosities in the 0.5-10 keV bin are re- 
ported against the B-band absolute magnitude (Mb) as 
open diamonds in Figure 1111 The comparison of the lu- 
minosities for the IR-detected (filled circles) and the un- 
detected objects (open diamonds) reveal some tendency of 
the former to display larger luminosities in both bands. Our 
high-z spheroidal galaxies tend to show larger values of the 
X-ray luminosity than the local counterparts in Figur jTTl 
(the local sample of spheroidal galaxies is taken from Ellis 
et al., 2006, small crosses in the Figure). In fact, the median 
of the X-ray luminosity for our high-z sources is on the order 
of log(Lx) ~ 42.2 erg/s, while the median of the X-ray lu- 
minosity for the local sample is log(Lx) ~ 40.5 erg/s. Even 
if the selection of our flux limited sample (mag(z) < 22.5) 
should ensure the completeness at optical wavelengths, to 
account for the possibility that our sample might be bi- 
ased towards intrinsically brighter systems (simply because 
fainter objects at high-z could fall out of the sample), we 
have also limited the comparison to the brighter systems. 
If we consider sources brighter than Mb\-21 (as marked by 
the vertical dashed line in Figure ITT]) . the median of the X- 
ray luminosity for our high-z sources is still on the order of 
log(Lx) ~ 42.2 erg/s, while the median of the X-ray lumi- 
nosity for the local sample is brighter than the previous case, 
log(Lx) ~ 41.3 erg/s. However, also at brighter optical lu- 
minosity the high-z; spheroidal population seem to populate 
brighter X-ray luminosities. However, the agreement with 
local data start to look better at even brighter optical lumi- 
nosities (Mb[-22). In this case, again log(Lx) ~ 42.2 erg/s 
for the high-z sample, but log(Lx) ~ 41.7 in the local Uni- 
verse. 

To derive the absolute B magnitudes we have computed 
the observed B magnitudes in the rest-frame, for each ob- 
ject, by convolving the corresponding best-fit redshifted- 
SEDs with the B-band filter. We have then converted the 
observed rest-frame magnitudes into their absolute values 
through the corresponding luminosity distances. Error bars 
in the Figure accounts for the uncertainties on the photome- 
try and on the 7\-correction applied to derive the rest-frame 
absolute B magnitudes. 

The X-ray luminosities have been proposed to allow 



Hidden Activity in High-redshift Spheroidal Galaxies 17 




Figure 11. Plot of the X-ray luminosity in the 0.5-10 keV bin versus the B-band rest-frame luminosity for the sample sources. Open 
diamonds are the objects with X-ray counterparts but no detectable flux at 24 (im, filled circles are X-ray sources with 24 (im detection. 
Source identification numbers refer to the ID column in Tabic 3 and are reported to the left and to the right of the data-points for 
diamonds and circles, respectively. The small symbols correspond to optical and X-ray luminosities of a local spheroidal galaxy sample 
by Ellis et al. (2006). 



further rough discrimination between AGN- and starburst- 
dominated phenomena: the value of Lx — 10 42 erg/s has 
been found to be a limit only rarely approached and never 
exceeded by star-forming galaxies (Franceschini et al. 2003; 
Persic et al. 2004). 

An important point about the origin of the X-ray en- 
ergy source is to understand how much of it might be related 
with the presence of long-lived emissions by hot plasma coro- 
nae or by old evolved stellar populations, instead of tracing 
ongoing activity phenomena. 

Discovered by the Einstein Observatory (Forman et 
al. 1979), the hot coronae around local early- type galaxies 
were interpreted as plasma distributions with masses up to 
several 10 10 Mq in quasi-static equilibrium in the galaxy 
dark matter halo. The heating and cooling mechanisms and 
timescales are not completely clear. The former, in partic- 
ular, might be related with residual star formation and su- 
pernovae explosions (e.g. Forman et al. 1994) . 

Because of these uncertainties, the interpretation of 
the X-ray emissions in our high-redshift spheroidal galax- 
ies might turn out to be ambiguous. However, the tempera- 
tures observed in local objects for such gas coronae are low, 
T ~ 10 7 K (Forman et al. 1994), consistent with the modest 



gravitational fields confining them. For the 8 galaxies de- 
tected in hard (2-8 keV) X-rays, our observed X-ray fluxes 
require much hotter plasmas (T > 10 8 K, computed by as- 
suming that only a single thermic component contributes to 
their X-ray emission), also taking into account the appre- 
ciable i\-correction. How much of the hard X-ray fluxes by 
early-type galaxies might be due to old stars in low-mass X- 
ray binaries (LMXB) is a problem already addressed in de- 
tail in a number of papers. Forman & Jones (1989), Fabbiano 
ct al. (1989) argue that only X-ray luminosities lower than 
3 x 10 40 erg/s can be contributed by LMXB's, while those 
higher than 10 41 erg/s are dominated by gaseous emission. 
For our typical sample sources, the X-ray luminosities ex- 
ceeding a few times 10 41 erg/s indicate that the LMXB con- 
tribution should be negligible, with the exception of source 
#6x, likely contributed by either long-lived plasma emission 
or LMXB's. 



For the other 5 objects in our sample undetected in the 
hard X-ray band the interpretation is less clear, and will be 
further discussed later. 
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Figure 12. Colour-colour plot with Shk j 54.5 versus Shk/Sr. Srk and Sr are the flux densities in the near-IR H K-hand and the 
optical R-band. The solid lines delimit regions of this plane typically populated by type-f AGNs, (lower region), type-2 AGNs (middle 
region), while galaxies are mainly distributed in the upper area of the plot (Franceschini et al. 2005; note that these latter use the 
ff-band flux instead of the present HK, hence a small shift in the lines by an average factor 1.3 has been applied to the delimiting lines 
here). In the upper panel the different symbols indicate the spectroscopic classification based on the optical spectrum. In the lower panel 
we differentiate spheroidal galaxies with 24 (im flux detection only (open squares), from those detected in both MIR and X-rays (filled 
squares), and from the X-ray emitting spheroids with no detectable MIR emission (crosses). For these latter the ID number is reported 
in the upper panel, while ID numbers in the lower panel refer to the sources detected in the X-rays but not in the MIR. The typical 
error bar is indicated in the upper left angle of the top panel. 



4.5 Colour-Colour Diagnostics and Broad-Band 
X-ray-to-IR Flux Ratios 

We have looked for additional colour diagnostics to further 
discriminate the nature of the emissions in our sample of 
spheroids. A combination of optical to near-IR colours has 
been recently suggested by Franceschini et al. (2005) as a 
diagnostic tool to disentangle among different AGN cate- 
gories and normal galaxies. The corresponding plot for our 
source sample is reported in Figure 1121 showing the ratio of 
the H K' -band over the 4.5 /im fluxes against the HK' to 
optical J?-band flux ratio. From bottom-left to top-right, the 
three domains delimited by the solid lines correspond to the 



three different populations: type-1 AGNs, type-2 AGNs and 
normal/starburst galaxies, in the order. The three regions in 
the figure have been empirically calibrated by Franceschini 
ct al. (2005) by using a sample of sources with known spec- 
tral types. Type-1 AGNs and optical quasars exhibit stan- 
dard blue optical colours (attributed to an accretion disk) 
and steeply rising near-IR spectra (due to the presence of 
dust re-radiation; e.g. Fritz et al. 2006), and are expected to 
populate the bottom-left corner of the plot. Normal or star- 
bursting galaxies have, according to the redshift, either blue 
optical and blue near-IR colours (at low-a), or red optical 
and red near-IR (at high-z), and tend to fall in the comple- 
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Figure 13. The ratio of the 4.5 /im band over the HK' fluxes against the ratio of the broad-band fluxes in the 0.5-8 keV band over the 
3.6 fim band. Different symbols correspond to different spectroscopic classification. The meaning of the symbols is the same as reported 
in Figure [12] 



mentary region of the plot outside the solid lines. Finally, 
type-2 and obscured AGNs populate intermediate colour re- 
gionfl 

In the upper panel the different symbols indicate the 
spectroscopic classification based on the optical spectrum. 
In the lower panel we differentiate spheroidal galaxies with 
24 /im flux detection only, from those detected in both MIR 
and X-rays, and from the X-ray emitting spheroids with no 
detectable MIR emission. 

As expected, our high- 2 spheroidal galaxies with either 
IR excess or X-ray emissions almost completely avoid the 
colour region corresponding to type-1 quasars (bottom-left 
corner). 



8 In their original version Franceschini et al. used standard John- 
son ff-band, instead of the HK 1 used here: solid lines in Figure 
[12] have then been corrected to account for this difference. 



Objects with both MIR and X-ray bright detections 
(filled squares in the bottom panel) dominate the type-2 
AGN regions. Top panel also reveals that the same region is 
occupied by sources with characteristic AGN lines (filled cir- 
cles). Sources lacking an X-ray counterpart (open squares) 
or those with an X-ray detection only (crosses) are mostly 
distributed in the region expected for normal/starforming 
galaxies. Sources detected at 24 /im and in X-rays fall mostly 
inside the boundaries of type-1 and type-2 AGNs. 

We observe a generally fair agreement of this diagnostics 
with that of Figure [9] Source #8, which falls in the AGN 1 
region in the plot of Figure 1121 is only marginally consistent 
with an AGN classification based on its hard-X flux and R- 
band magnitude. However, this object is spectroscopically 
classified as ELG, implying a difficult interpretation of its 
MIR excess. 

The near-IR colour information is complemented in Fig- 
ure[l3]with that on the broad-band X-ray-to-IR flux ratios, 



20 G. Rodighiero et al. 



following a diagnostic scheme also discussed in Franceschini 
et al. (2005). The two panels report data with different sym- 
bols corresponding to different classification criteria for the 
sources, as detailed in the inserts. The bottom panel includes 
also sources with X-ray detection only (no MIR). Bona-Sde 
AGNs, with clear excess fluxes in X-rays and in the MIR (4.5 
/im), occupy the region indicated in the top-right quadrant, 
while objects with less extreme X-ray and IR emissivities, 
hence more likely dominated by stellar activity, are spread 
in the other quadrants of the plot. 

We see an overall good level of agreement among these 
various classification tools, with the exception of those based 
on the optical emission/absorption lines providing a classi- 
fication often at variance with the others. 



5 DISCUSSION 

5.1 Local versus High-Redshift Spheroidal 
Galaxies 

With few exceptions, early-morphological type galaxies in 
the local universe appear as typically inactive systems, al- 
most completely devoid of an interstellar medium, with no 
relevant star-formation or AGN activity in place. This view 
has been essentially confirmed by a variety of observations, 
including those carried out with space IR and X-ray ob- 
servatories. The IRAS survey first revealed that a fraction 
of local spheroidal galaxies contain some amounts of ISM 
gas and dust (Jura et al. 1987; Knapp et al. 1989). The In- 
frared Space Observatory (ISO), with better sensitivity and 
higher resolution, has shown the presence of trace amounts 
of dust and Polycyclic Aromatic Hydrocarbon molecules 
(PAHs) in the MIR spectra of a few local early-type galaxies 
(i.e. Madden et al. 1999; Athey et al. 2002; Xilouris et al. 
2004), results also recently confirmed by observations with 
the Spitzer Space Telescope (Pahre et al. 2004). Spitzer IRS 
spectra reported by Bressan et al. (2006) show that, in the 
bulk of local early- type galaxies, the MIR emission is due to 
circumstellar dust shells produced by AGB stars. 

In such local populations the dust contribution to the 
MIR light is quite modest, however, providing a tiny frac- 
tion of the bolometric galaxy emission. The corresponding 
amount of dust is negligible. 

On the contrary, we have found that a significant frac- 
tion (~15%) of high-redshift spheroidal galaxies are detected 
with substantial mid- and far-IR luminosities, which are a 
factor ~10 higher than found in the local counterparts (see 
Figure [6} . The large MIR excess emissions are apparent in 
the SED's of Figs. [4] and [5] particularly evident if compared 
with the integrated spectra of local objects reported e.g. 
by Athey et al. (2002) and Bressan et al. (2006). These 
emissions cannot be interpreted in terms of simple stellar 
mass loss. The latter interpretation of our observed fluxes 
would require both mass loss rates quite larger than ex- 
pected during any evolutionary phases of spheroidal galax- 
ies, and metallicities of the main stellar populations factors 
2-3 higher than solar (Bressan et al. 1998, 2001). Further 
evidence that circum-stellar dust by evolved stars cannot be 
responsible for the excess IR emissions is provided by the 
emission lines detected in all the optical spectra of figure [7] 
and by the large X-ray luminosities of the majority of the 
sample objects. 



A related question is about the origin of the observed 
X-ray emission for the sub-sample objects with no MIR de- 
tection. Less information is available for these spheroidal 
galaxies, and correspondingly more difficult is their interpre- 
tation. These sources tend to be less active than the MIR 
emitting ones, as for example indicated by their typical X- 
ray luminosities (log(Lx) < 42 erg/s), and by the fact that 
5 out of 10 of them do not show appreciable emission lines 
in Figure [5] Here again, however, the properties for the ma- 
jority of the sources appear to require the presence of an 
activity source, as summarized in Table 1. 

5.2 Nature of the IR and X-ray Emissions in 
High-z Spheroidal Galaxies 

5.2.1 The IR-loud Population 

With reference to the 19 spheroids with detectable MIR 
emission, we have found 8 objects for which all our diag- 
nostic diagrams provide consistent evidence for a dominant 
AGN presence. These are sources #0, 2, 3, 4, 7, 9, 13 and 
20. Only four of them (#2, 4, 7, 13) display type-2 AGN 
features in the optical spectrum. 

Source #8 appears as a transition object in all diagrams 
and the SED seems to require a dusty-torus contribution to 
the MIR. Similarly, source #6 shows a LINER spectrum and 
its SED suggests an AGN contribution in the MIR. Source 
#10 also appears as a transition object, with colours in Fig- 
ure [T2] close to those of type-2 AGNs and a large X-ray 
luminosity in Figure [TT1 

For eight further sources (#1, 5, 11, 12, 15, 16, 17, 19), 
none of our diagnostic diagrams indicates the presence of 
an AGN. For most of these, instead, the rest-frame ~12 
fim and X-ray luminosities in Figs. [6] and [11] require the 
presence of star formation activity in the galaxy. Based on 
our spectral best-fits to the SED data, we have computed the 
bolometric luminosities and from these estimated the star- 
formation rates (SFR's) using recipes discussed in Rowan- 
Robinson et al. (1997) and Elbaz et al. (2002). We infer 
values of SFR from ~ 1 to ~ 20 M Q /j/r. For the two low- 
redshift galaxies #5 and #15, the SFR is instead around 0.2 
M & /yr at most. 

5.2.2 The Population Emitting in X-rays only 

A basic uncertainty is whether an activity (AGN or star- 
burst) is present at all in these galaxies, or if their X-ray 
flux might be entirely explained by long-lived emissions by 
hot plasma coronae or LMXB populations. Eight out of the 
14 X-ray-only galaxies are detected in the hard X-ray (2- 
8 keV) band with luminosities larger than 10 41 erg/s. As 
discussed in Sect. 14.41 the combination of the large X-ray 
luminosity and hard X-ray band implies either an AGN or 
an ongoing starburst to reside in these sources. Of these, 
at least 5 (ID #3x, 5x, 9x, llx, 12x) show evidence for an 
AGN in the large value of the X-ray luminosity (Figure [T3)l . 
the large X/O ratios (Figure O, and in the optical spec- 
tra (Figure [HJ. Sources #4x and #14x are likely to host a 
luminous starburst, while source #7x is a low-optical lumi- 
nosity, moderately active galaxy in X-rays, and could host 
a low-luminosity AGN. 

Two of the remaining 6 objects detected only in the soft 
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X-rays (ID #lx and #13x), show prominent [Oil] and other 
emission lines and are star forming galaxies. Source #2x is 
a very luminous soft X-ray and very bright galaxy in the 
optical (Figure II ip . with a passive (absorption-line, post- 
starburst) optical spectrum, and likely to host an obscured 
starburst, or post-starburst young stellar population. 

Source #6x is, on the contrary, a very low-luminosity 
local object in both soft X-rays and the optical. The X-ray 
emission is likely to originate from a hot plasma corona, also 
confirmed by the complete absence of optical emission lines. 

Finally, for the soft X-ray emitting objects #8x and 
#10x we lack an optical spectrum. Both are moderately lu- 
minous in the optical and X-rays (Mb — —22, Lx — 5 x 10 41 
erg/s) and all their properties are consistent with those of 
starbursts. 

5.3 Galaxy Morphologies 

The unique multi-band ACS imaging add valuable physical 
information on the sample galaxies. Indeed, if the activity 
is triggered by merging or interactions, we expect we would 
be able to detect it. 

The images in Figs. [T] and [2] first confirm that the far-IR 
emitting objects appear to have been detected in a higher 
activity stage than the X-ray emitting only. Of the 19 IR- 
loud galaxies, 12 show morphological peculiarities or some 
evidence of interactions, while only 4 of the 14 sources de- 
tected in X-rays show morphological irregularities. 

In some cases (ID #4, #8, #9, #10, #11, #9x) the 
spheroidal galaxies look part of a group. 

However, the majority of the X-ray emitting and a good 
fraction of the IR-loud objects (i.e. 15 out of 33 in the total 
sample) appear well isolated and showing the regular ellip- 
tical contours of spheroidal galaxies. 

5.4 Statistical Considerations 

From the combined diagnostics considered in this work 
(morphology, spectroscopy, X-ray emissions, colours, broad- 
band luminosity ratios), a rather heterogeneous scenario 
emerges concerning the nature of the population of "active" 
spheroidal galaxies at redshifts from 0.3 to 1.2. 

In summary, among the 33 spheroids with excess MIR 
or X-ray emissions, 13 show clear evidences for the presence 
of an obscured AGN as mainly responsible. For 16 objects (8 
among the IR-loud and 8 more in the X-ray loud classes) the 
most likely interpretation is in terms of an hidden starburst, 
with typical SFR values between 1 and few tens M^/yr. The 
remaining 4 galaxies display less characteristic, intermediate 
physical properties. Among them, the rate of star-formation 
appears to be modest in the 2 lowest-redshift galaxies (#5 
and 15) or even absent in #6x. 

In conclusion, almost 20% of the optically-selected pop- 
ulation of faint high-z spheroidal galaxies show evidence of 
activity at the instant of the observation. Roughly half of 
them (~10% of the total optically-selected spheroidal galaxy 
population), appear dominated by AGN, the other half by 
starburst activity (another ~10% of the total). In either 
case, the duration of the phenomenon is expected to be 
short, of the order of 10 8 yr (e.g. Forster-Schreiber et al. 
2003; Cavaliere & Padovani 1989), and likely to follow merg- 
ing or interaction with a gas- rich system, able to funnel fresh 



gas into the galaxy nucleus. The process would then be ex- 
pected to contribute new stars and to increase the nuclear 
super-massive black hole in the AGN sub-class. 

Compared with the 5-7 Gyrs corresponding to the Hub- 
ble times at the source redshifts, the expected short duration 
of the processes and the significant fractions of the objects 
detected in the active phase are consistent with the fact that 
the whole spheroidal population undergoes few to several 
such episodes at z ^ 1. 

Then assuming for example as a reference for the SFR 
the value of 10 M^/yr and a star-formation timescale of 
10 s j/r, this would imply several such episodes for each 
galaxy, each one contributing of the order of few percent 
of the final stellar content of the galaxy on average. For the 
most massive galaxies this is not expected to modify sig- 
nificantly the stellar mass below z ~ 1, while for objects 
with M < 10 11 Mq the increase with cosmic time would be 
expected to be more substantial. 

Altogether, these results are at least not inconsistent 
with recent analyses of the evolution of the galaxy stellar 
mass as a function of redshift based on complete near-IR se- 
lected samples with spectroscopic follow-up (Dickinson et al. 
2003; Bundy et al. 2005; Fontana et al. 2004; Franceschini et 
al. 2006). Moreover, starting from a similar (and partially 
overlapping) morphological sample as used in this paper, 
Treu et al. (2005) confirmed the presence of star formation 
in spheroids at high redshift. They studied the spectroscopic 
properties of spheroids and the evolution of the fundamen- 
tal plane, and confirm a scenario in which the efficiency of 
SF is enhanced in the most massive galaxies at the higher 
redshifts, and proceeds with longer timescales in the smaller 
ones (the downsizing effect). 

Our results, in which AGN and starburst activities are 
tightly intertwined and actually almost inextricable in sev- 
eral of our galaxies, emphasize the continuity of the two 
processes and further support the case for a co-evolution of 
quasars and bulges (Franceschini et al. 1999; Granato et al. 
2004). 



6 CONCLUSIONS 

To complement searches for recent or ongoing star-formation 
and nuclear accretion in high-redshift spheroidal galaxies, we 
have exploited deep MIR and X-ray observations by Spitzer 
and Chandra in the GOODS-N area. This allowed us to 
identify signs of hidden activity that might escape optical 
investigations based on colours or strong morphological sig- 
natures of merging or interactions. 

Our reference sample are 168 morphologically classified 
spheroidal (elliptical/lenticular) galaxies with zab < 22.5 
selected from GOODS ACS imaging. We have found that a 
significant fraction of the E/SO's show unexpected emissions 
in the MIR: nineteen of them have 24-pm detections in the 
GOODS catalogue, while 6 have detectable radio flux at 1.4 
GHz with >40 fiJy. 

We have also looked for traces of hidden activity in 
our optically-selected spheroidal population in the 2 Ms 
Chandra catalogue and detected significant flux from 25 
sources (13 of which having 2-8 keV hard X-ray emission). 
Of these, 11 belong to the MIPS 24-^im catalogue and 14 
more are detected in X-rays only. 
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The nature of the near- and MIR emission observed by 
Spitzer is analyzed through the SEDs based on the available 
multi-wavelength photometry, including X-ray, UV, optical, 
near-, MIR, and radio fluxes. The inferred amount of dif- 
fuse dust appears substantially in excess of that expected 
by mass loss from late-type stars. We have used in addi- 
tion a variety of diagnostics, including optical spectroscopy, 
colour-colour diagrams, and broad-band flux ratios to inter- 
pret the nature of the energy source. 

When the available independent diagnostics are com- 
pared, in general they provide consistent classifications 
about the nature of the activity in the spheroidal population. 
However, the SED fitting study alone tends to privilege an 
AGN predominance in almost all sources. X-ray diagnostics 
(luminosities and spectral analysis) and colour-colour dia- 
grams instead provide much correlated classifications. Given 
that, in principle, none of these diagnostics alone can be 
considered as conclusive, only trough a panchromatic com- 
parison of them we can reach an accurate comprehension of 
the underlying physical processes. 

Our multi-wavelength analysis of the X-ray and MIR 
properties leads us to conclude that at least 8 of the 19 
IR-detected sources should hide an obscured AGN, while 
the X-ray undetected sources are more likely dominated by 
ongoing star formation. 

More uncertain is the interpretation of spheroidal galax- 
ies detected in the X-ray bands only and undetected in the 
MIR, a basic question being to understand whether an ac- 
tivity source is present at all, or if the X-ray fluxes might 
be entirely explained by long-lived emissions by hot plasma 
coronae or LMXB populations. Our analysis shows evidence 
for the presence of an AGN or an ongoing starburst in the 
majority of them. 

We conclude that 15% to 20% of the original sample 
of 168 spheroidal galaxies in GOODS-N are detected during 
phases of prominent activity, probably following an event 
of merging/interaction. Our morphological analysis has re- 
vealed that only roughly a half of the galaxies show evi- 
dence of some peculiarities, the other half displaying stan- 
dard spheroidal morphologies. This might be taken as evi- 
dence that our IR/X-ray selection tends to pick up the latest 
stages of a merging/interaction process. 

Thirteen galaxies (8% of the complete sample) show 
evidence for the presence of an obscured AGN, and for 
other 16 (10%) the most likely interpretation is in terms 
of a hidden starburst. Due to the expected short lifetimes of 
these IR and X-ray emissions, this observed fraction implies 
widespread activity in this class of galaxies during the cos- 
mic epochs - z ~ 0.3 to z ~ 1 - subject to our investigation. 
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* ACS data from Bundy et al. (2005) 



